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Braulio V. Sanchez 
ABSTRACT 

An investigation has been conducted to ascertain the response of the rotational motion of the 
Earth to forcing functions produced by tie water mass redistribution due to the ocean tides. In 
particular, the components of displacement of the rotation axis at the surface of the Earth were 
obtained. The investigation also addressed the larger question concerning the possibility of excita- 
tion of the Chandler wobble of the Earth. 

In general, the results show the existence of a polar wobble as a response to each of the com- 
ponents of the ocean tides. The magnitude of the polar displacement depends on two factors; the 
amplitude of the tidal component and its period (in relation to the Chandler period). 

The maximum periodic contributions are: the Doodson’s component number 055.565 with a 
period of 18.613 years and 50 cm of polar displacement, the annual component 056.544 with 37 cm 
of polar displacement and the semi-annual 057.555 with 32 cm. The tidal components with daily 
and semi-daily periods yield very small polar displacements of the order of 0.01 cm. The combined 
effect of all the periodic components can yield as much as 90 cm of pole displacements. 

The changes produced by the ocean tides in the products of inertia are periodic and regular, 
therefore, they cannot be the source of excitation of the Chandler wobble. 
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OCEAN TIDAL EXCITATION OF POLAR MOTION 


INTRODUCTION 

The objective of this investigation is to ascertain the response of the rotational motion of 
the Earth to forcing functions produced by the water mass redistribution due to the lunisolar 
ocean tides. In particular, it is desired to obtain the components of displacement of the rotation 
axis at the surface of the Earth. The investigation also addresses the larger question concerning 
the possibility of excitation of the Chandler wobble of the Earth. The problem concerning the 
interaction between the angular velocity of rotation of the Earth and the induced oceanic tidal 
response has been treated by others, i.e.: Haubrich and Munk (1959), Dahlen (1976), and it will 
not be dealt with here. 

One of the areas of interest in present day geophysics deals with the excitation of the polar 
motion of the Earth, that is, how do the different physical phenomena taking place inside and on 
the surface of the Earth affect the position of the rotation axis and the rate of rotation. 

Among the many physical phenomena one must consider the redistribution of the mass of 
the ocean produced by the gravitational tides due to the Sun and the Moon. Such a redistribu- 
tion of mass yields changes in the inertia tensor of the Earth which themselves are the causes of 
changes in the position of the rotation axis and the rate of rotation. All the changes involved 
are functions of time with various periodicities. 

In order to ascertain the ocean tidal effects the Liouville equations of motion for a nonrigid 
body have been simplified by neglecting higher order terms and an analytic solution has been 
found for the cUjj and Wy components of angular velocity. Such a solution involves the magni- 
tude of the changes in the products of inertia I^z and ly^ due to the various ocean tidal compo- 
nents as well as the frequencies of the tides and the natural frequency of rotation of the Chand- 
ler wobble. 


1 



The evaluation of the changes in I^z and lyz can be divided in two categories. Changes pro- 
duced by ocean tides with periods of the order of a day and changes produced by long period 
ocean tides. The method of evaluation depends on the category and the effects produced by 
each are considerably different in order of magnitude. 

The products of inertia due to the short period ocean tides have been obtained by making 
use of spherical harmonics representations of the tide heights. These representations are fits to 
numerical solutions of the Laplace tidal equations and as such represent non-equilibrium solu- 
tions. In particular, the harmonic coefficients used were those provided by C. Goad 0978), the 
LTE solution used is due to R. Estes (1977). 

The long period ocean tides were treated differently since no harmonic expansions were 
available. The equilibrium response of an ocean covering the entire surface of the Earth was 
formulated, the tidal height for the real Earth is then obtained by means of the ocean function 
which is equal to one over the oceans and zero over the continents. In particular the (8 x 8) 
ocean function harmonic expansion due to Munk and MacDonald (1975) was adopted. In order 
to solve for the products of inertia it is then necessary to solve the integrals over the sphere of 
the product of three surface harmonics. Such a solution can be expressed in terms of the “3 - j” 
symbols often found in quantum mechanics. 

In order to represent the lunisolar tidal potential the expansion by A. T. Doodson (1954) 

was adopted with the extensions and corrections due to Cartwright and Taylor (1971) and Cart- 

% 

Wright and Edden (1973). 

Finally, the solutions for the and Wy components of angular velocity were mapped into 
components of displacement of the rotation axis at the surface of the Earth. 

1. SOLUTION TO THE LIOUVILLE EQUATIONS 

The Liouville equations of motion were first given by Liouville (1858). For purposes of 
this investigation the following assumptions can be made: 
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(1) the external moments are zero, 

(2) the relative angular momentum terms vanish, 

(3) the moments of inertia are constant and considerably larger than the products of inertia, 

(4) the equatorial moments of inertia are equal, 

(5) the component of angular velocity is a constant and much larger than and cOy. 


Then, neglecting products of small quantities the following equations result. 

Au}]( “ (A - C) cOj Wy — <*^2 ~ ^yz ^z 

AtOy + (A - C) COj, - iyz ^xz ^2 
Equa* ons ( 1 . 1 ) can be written as follows 

+ n^cux = A-‘ { I'xz lOz - nco^ - iy^ + nwz)} 
o3y + n^Wy = A-* (ly 2 ~ nco^ (<j2 + ncoz)} 

where 

n = (C - A)w 2/A 
Let 

^XZ ~ ^XZ “ ^xz) 

*yz ^ ^yz (ft - 0y2> 

Then 

a>x = K, cos (ft - 0xz) + ^2 sin (ft - ^y^) 
cOy = K3 cos (ft - 0y2> + K4 sin (ft - 0^2) 
iCi = -M^z cjz(f2 + nwz)/[A(n2 - f2)] 

K2 = Myz <Oz f(cjz + n)/[A(n 2 - f 2 )] 

K3 = (Myz/Mxz)K, 

K4 = - (M„/Myz)K2 


(1.1) 


( 1 . 2 ) 


(1.3) 


(1.4) 


(1.5) 


Equations (l.S) can be written as follows, 
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Wx = Wx cos (ft - 
Wy “ Wy cos (ft “ 0<jy) 

Wx= {Kj + +2K, Kj 

Wy = {Kl + Kj + 2 K 3 K 4 

^tox “ {(^1 ^xz 

= arc tan ((K 3 sin ^y^ 

2. TIDAL POTENTIAL 

The tidal potential can be written in different ways. In particular, the second degree zonal 
component in spherical coordinates is given by 

v 5(R) = (R)2 P5(cos0)q; (2.1) 

R: radius of the Earth, 

P® : second degree legendre polynomial, 

; a function which depends on the disturbing body. 

Doodson's expansion for the long period lunisolar potential is given by (Bretreger, 1978), 

V°(R) = -Gp pOfcosd) ^ Aj cos(Oj) (2.2) 

i 

Comparing Equations (2.1) and (2.2) it follows that 

q° = -(Gd/R 2) Aj cos(Oi) (2.3) 

i 

where Gq is Doodson’s constant and Aj, Oj are the amplitudes and arguments for the various 
waves (see Appendix 2). 

The principal terms of the low frequency tides are given by Cartwright and Edden (1973). 
For the purposes of this investigation the following terms will be considered: The 17 waves of 


sin ( 0 XZ ” ^yz^j 

, 1/2 

sin (0y2 - ^xz)} 

+ Kj cos 0yj,)/(Kj cos ~ *^2 ^yz)} 

+ Iv^ COS *^yz ” ^4 ^xz^} 
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group (0,5) with periods ranging from 91 days to infinity. The 2 terms with largest amplitude in 
group (0,6) with periods of 27 and 31 days. The term with largest amplitude in group (0,7) with 
period of 1 3 days. The group terminology is that of Doodson. 


Table (2.1) below gives the tidal components with their amplitudes and periods. 

Table (2.1) 

Tidal Components 


Doodson’s No. 

Period (days) 

Amplitude (Aj) 

055.555 

OO 

0.73806 

055.565 

6798.2530 

-0.06556 

055.575 

3399.1265 

0.00064 

055.765 

1305.4098 

-0.00009 

056.544 

385.99876 

0.00009 

056.554 

365.25964 

0.01156 

056.536 

365.22476 

-0.00062 

056.564 

346.63545 

-0.0001 1 

057.345 

212.32434 

-0.00005 

057.355 

205.89383 

0.00074 

057.555 

182.62110 

0.07281 

057.553 

182.62982 

0.00029 

057.565 

177.84369 

-0.00180 

057.575 

173.30987 

-0.00040 

058.554 

121.74934 

0.00426 

058.564 

119.60730 

-0.00007 

059.553 

91.312730 

0.00017 

063.655 

31.811921 

0.01579 

065.455 

27.554564 

0.08254 

075.555 

13.660791 

0.15647 
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3. EQUILIBRIUM TIDES 

The surface of an ocean covering the entire surface of the Earth is raised with respect to the 
ocean bottom by the following amount: 


(1 -f k - h) 
g 


V'*' 

Vj, 


k, h: love numbers, 

g: gravitational acceleration, 


where V| denotes the disturbing potential evaluated at the surface of the Earth. 


(3.1) 


Let n denote the surface area covered by the oceans and define as follows; 

k - h) 




g 


Vjds, 


(3.2) 


where the surface integral is taken over the area of the continents. represents the quantity 
which must be added to {q in order to satisfy conservation of mass. Note that, 



(1 + k - h) 
g 


V] ds = 



(1 + k - h) 
g 

(1 + k - h) 
g 


V] ds 


V] ds 


(3.3) 


The first term on the right hand side of Equation (3.3) goes to zero, therefore. Equation (3.2) 
yields: 





//. 


Oceans 


(1 -t-k -h) 
f 


VT ds. 


(3.4) 


The resulting height of the ocean surface is then given by; 

{ = to + (3.5) 

Inserting Equations (3.1) and (3.4) into Equation (3.5) yields; 
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(1 ^-k-h) 
g 




Vj ds 


Oceans 


Let the “ocean function” be defined by: 


°° n / air cos m ij/ ' 

= y Y pir(cosd) " 

^ sin mi//, 

n*0 m*0 ' n 


and the distu'bing potential by: 


Vj = R2 ^ Pin{cos0) 


q™ cos m 


sin m ^ , 


m»0 


(3.6) 


(3.7) 


(3.8) 


Furthermore, consider only the terms up to n = 9 in Equatio»i (3.7) and m = 0 in Equation (3.8). 
Then, 


n 


- I f ds = 4ff a”, 

•' Sphere 


// fj 

•'Oceans •' *''5 


Sphere 


f(0,i//) ds = ^ q®. 


(3.9) 


(3.10) 


Equation (3.6) can then be written as follows: 




P»(cos9) - ^ 


a§ 




The contribution to the products of inertia l^y and lyj, are given by: 


‘XZ 


*yz 


■//. 

■// 


f(0,\//) XZ dm. 


Sphere 


f(0,^{/) yz dm. 


Sphere 


where 


X = R sin0 cos^ 
y = R sin0 sini/» 
z = R COS0 

dm = (R d0 • R sin0 di//) . 

Py, : density of sea water 


(3.1 n 


(3.12) 


(3.13) 


(3.14) 
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Making use of Equations (3.14), (3.) 1) and (3.7) allows Equations (3.12) and (3.13) to be written 
as follows; 


‘xz 


.(1+k-h)/'^’^/’’^’'' _ /aj" cos m vl/ \ 

= r6 L / y y p|;\cos0) " 

8 Jq Jq ^ Z-rf \bjr sin m / 

^ " n=0 m»0 ' ” 


a® 1 1 

P^(cos0) — ^ q® •— P^(cosO) sin0 cos\l/ d0 d^, 
Sap 3 


lyz " Pvf R* 


(1 + k - h) 


'I I?, U .n . 


g 


(3.15) 


• r 


P9(cos0) - 


5a® 


q9 •— P*(cos0) sin0 sini// d0 di^^. 


Integration of Equations (3.15) yields the following (see Appendix 1): 


^xz ~ ^ ^xz ^l2 
lyz = K Ty, q® 


(3.16) 


K = (47r/g) + k - h) 

T„ = (l/35)a; + (2/21)ai - (l/25)(a® aj/aj) 
Ty^ = (l/35)b; + (2/21)b{ - (l/25)(a; bj/a®) 

Making use of Equation (2.3) yields 


*xz “ Txz Ou R*^ ^ Aj cos(o£j) 

i 

Iy2 = -K Tyj Gp R'^ ^ Aj cos(ai) 

i 


(3.17) 


4. NON-EQUILIBRIUM TIDES 

The method of solution employed for the equilibrium tides becomes less and less applicable 
as the period of the forcing function becomes shorter. Specifically, the equilibrium theory is 
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apprc nate to the extent that (P^/Pf )2 <5C 1 , where P„ is a natural period of the system and Pj- 
is the period of the forcing function. 


The products of inertia due to the short period ocean tides have been obtained by means of 
existing spherical harmonics representations of the tide heights. These expansions are fits to nu- 
merical solutions of the Laplace tidal equations. 

The spherical harmonic coefficients for the tide heights were obtained from C. Goad (1978), 
the LTE solution used is due to R. Estes (1977). 

Table (4.1) gives the values of the coefficients used in this investigation. 


The expression for the tide height is given by 


^ P;f(cos0) 
n m=0 


/B!P cos fi t\ 

+ 1 I sin m \£/ 

\K sin ty 


r/Ajl'cosfi t’ 

L\Cn sin fi 

>} 


cos m\I/ 


The products of inertia are given by 


*xz 


lyz ~ 


= //. 


xz dm 


yz dm 


Therefore, making use of Equation (4.1), 

Ixz = (47 t/5)R^ Pw(^ 2 ^°sfi t + C* sinfj t) 
lyz = (4 it/5)R'* Pw(® 2 t + Dj sinfi t) 


(4.1) 


(4.2) 


(4.3) 


or 
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^xz ^xz “ *^xz) 

lyz ~ ^yz “ 0yz^ 

Mxz = a[(Ai)2 + (C^)2]i/2 

Myz = al(B\y^ + (D1)2]V2 (4.4) 

<}>xz ~ 3 fc tan (C2/A2) 

^y2 = arc tan (D^/Bj) 
a = (47t/5) R'* 


Table (4.1) 

Coefficients for Short Period Tides 


Doodson No. 

Darwin 

Period (hr) 

Ai(102) 

cj(io^) 

?4(102) 

dJ(102) 

273.555 

S2 

12.00 

-0.5642 

-0.2811 

-0.06755 

-0.3763 

255.555 

M2 

12.42 

-1 .274 

-0.8321 

-0.4521 

-0.9790 

245.655 

N2 

12.66 

-0.2296 

-0.1622 

0.1923 

-0.04628 

145.555 

Oi 

25.82 

0.3666 

3.070 

1.974 

2.322 

165.555 


23.93 

-0.7530 

2.110 

3.676 

3.961 


5. NUMERICAL RESULTS 

The results obtained by n^eans of Equations (3.17) and (4.4) are given in Table (5.1) on the 
following page. The magnitudes M^z, Myj. are given in units of 1 0^3 gnJ-cm^. The values used 
for the constants are the following: 

Gd = 2.627723 (10^) cm2/sec2 
R = 6.378388 (10«) cm 
p^ = 1.030 gm/cm3 
A = 8.016604490270 (10‘W)gm-cm2 
g - 980 cm/sec^ 

1 + k - h = 0.7 
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Table (5.1) 

Tidal Contributions to the deducts of Inertia, Epoch: 
1899, Dec. 31, 12h 


Doodson No. 

M,2 (10-») 

0XZ (deg) 

My, (10-33) 

0yz (deg) 

055.555 

-1500.00 

0 

-735.00 

0 

055.565 

133.00 

259 

65.30 

259 

055.575 

-1.30 

158 

-0.63 

158 

055.765 

0.18 

207 

0.08 

207 

056.544 

-0.18 

-260 

-0.08 

-260 

056.554 

-23.40 

-1 

-11.50 

-1 

056.556 

1.26 

200 

0.61 

200 

056.564 

0.22 

257 

0.10 

257 

057.345 

0.10 

-368 

0.04 

-368 

057.355 

-1.50 

-109 

-0.73 

-109 

057.555 

-147.00 

199 

-72.50 

199 

057.553 

-0.58 

-3 

-0.28 

-3 

057.565 

3.65 

98 

1.79 

98 

057.575 

0.81 

-3 

0.39 

-3 

058.554 

-8.65 

197 

^.24 

197 

058.564 

0.14 

97 

0.06 

97 

059.553 

-0.34 

196 

-0.16 

196 

063.655 

-32.00 

45 

-15.70 

45 

065.455 

-167.00 

-63 

-82.20 

-63 

075.555 

-318.00 

180 

-155.00 

180 

165.555 

9.59 

109 

23.10 

47 

145.555 

13.20 

83 

13.00 

49 

245.655 

1.20 

-144 

0.84 

-13 

255.555 

6.51 

-146 

4.62 

-114 

273,555 

2.70 

-153 

1.63 

-100 
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The values for the coefficients of the ocean function are those given by Munk and MacDon- 
ald (1975). The normalization factors were taken into consideration. 

The values for the amplitudes Aj are given in Table (2.1) for the equilibrium tides, for the 
short period tides Table (4.1) gives the values for the coefficients. 

The phases ^’yz correspond to 1899 December 31 at I2h Om Osec ephemeris time 

(see Appendix 2). 

Table (5.2) presents the tidal heights for the equilibrium tides as obtained from Equation 
(3.11) with 9=0° and cos («j) = 1. These values therefore represent an upper bound. 


Table (5.2) 

Maximum Equilibrium Tide Heights 


Doodson No. 

Tide Height (cm) 

055.555 

-14.30 

055.565 

1.27 

055.575 

-0.12 (lO'i) 

055.765 

0.17 (10-2) 

056.544 

-0.17 (10-2) 

056.554 

-0.22 

056.556 

0.12(10-1) 

056.564 

0.21 (10-2) 

057.345 

0.97 (10-3) 

057.355 

-0.14(10-1) 

057.555 

-1.41 

057.553 

-0.56(10-2) 

057.565 

0.35 (10-1) 

057.575 

0.77(10-2) 

058.554 

-0.82(10-1) 

058.564 

0.13 (10-2) 

059.553 

-0.33(10-2) 

063.655 

-0.30 

065.455 

-1.60 

075.555 

-3.04 
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The maximum x and y components of displacement for the position of the pole are given in 
Table (5.3). These values are obtained by means of the following relations, 

X = (Wx/w2)R 

y = (Wy/W2>R 

where W^, Wy are given by Equations (1.6) and 

= (2jr/86400) rad/sec. 

The value of “n” (Eq. 1.3) has been chosen so as to obtain a Chandler period of 428 days. 
These values represent maximum values since it has been assumed that cos (oj) = 1. Table (5.4) 
shows the variation in d = vx^ ^ function of “n” for the epoch Dec. 31 , 1899, 12h 

ephemeris time. The value of “n” has been chosen so as to obtain Chandler periods of 400, 428, 
444 and 460 days. 

The last row of Table (5.4) gives the total pole displacement for the sum of all the periodic 
components, that is the component 055.555 is not taken into account since its effect can be nul- 
lified by a proper choice of coordinates. 

Figures 5.1, 5.2 and 5.3 show the pole displacement due to the components 055.565, 
056.554 and 057.555 respectively, the epoch is Dec. 31, 1899, 12h. In each case the time span 
corresponds to the period of the tidal component: 6798 days, 365 days and 182 days 
respectively. 

Figures, 5.4 and 5.5 show the x and y components of pole displacement generated by the 
sum of all the periodic components (055.555 excluded), the epoch is the same as above, time 
steps of 67.98 days were used. 

Figure 5.6 is a plot of + y2 with x and y obtained as explained above. A Chandler 

period of 428 days has been assumed to generate the results shown in Figures 5.1 to 5.6. 
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Table (5.3) 

Maximum Components of Pole Displacement 
Chandler Period = 428 Days 


Doodson No. 

X (cm) 

y (cm) 

y/ + y2 (cm) 

055.555 

510.000 

250.000 

568.0000 

055.565 

45.500 

22.500 

50.7000 

055.575 

0.450 

0.220 

0.5000 

055.765 

0.070 

0.040 

0.0800 

056.544 

0.300 

0.320 

0.4300 

056.554 

24.800 

27.300 

36.8000 

056.556 

1.320 

1.460 

1.9600 

056.564 

0.170 

0.190 

0.2500 

057.345 

0.010 

0.020 

0.0220 

057.355 

0.220 

0.330 

0.3900 

057.555 

17.200 

26.900 

31.9000 

057.553 

0.060 

0.100 

0.1100 

057.565 

0.440 

0.690 

0.8100 

057.575 

0.100 

0.160 

0.1800 

058.554 

0.520 

0.920 

1.0500 

058.564 

0.008 

0.010 

0.0120 

059.553 

0.010 

0.020 

0.0220 

063.655 

0.410 

0.810 

0.9000 

065.455 

1.850 

3.700 

4.1300 

075.555 

1.730 

3.470 

3.8700 

165.555 

0.010 

0.010 

0.0140 

145.555 

0.010 

0.010 

0.0140 

245.655 

0.001 

0.001 

0.0014 

255.555 

0.006 

0.005 

0.0078 

273.555 

0.002 

0.002 

0.0028 
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Table (5.4) 



Pole Displacement vs. Chandler Period, 
Epoch: 1899, Dec. 31, 12h 



v/ x2 + y2 (cm) 


Doodson No. 

400 Days 

428 Days 

444 Days 

460 Days 

055.555 

531.0000 

568.0000 

590.0000 

611.0000 

055.565 

9.3000 

10.0000 

10.4000 

10.8000 

055.575 

0.4300 

0.4600 

0.4800 

0.5000 

055.765 

0.0640 

0.0690 

0.0730 

0.0760 

056.544 

0.9100 

0.3300 

0.2500 

0.2100 

056.554 

41 .9000 

23.9000 

19.4000 

16.4000 

056.556 

2.2700 

1.3100 

1 .0700 

0.9100 

056.564 

0.2 ’00 

0.1900 

0.1700 

0.1500 

057.345 

0.0140 

0.0130 

0.0120 

0.0110 

057.355 

0.3500 

0.3400 
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Figure 5.1. Pole Displacement, Component 055.565 
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Figure 5.2. Pole Displacement, Component 056.554 
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Figure 5.6. Total s/\^ + vs. Time 





6. CONCLUSIONS 


The validity of most of the results is predicated upon the reality of the equilibrium response 
of the oceans to the long period components of the lunisolar potential. Within the context of 
this assumption and the simplifications made, the results indicate the following. 

There exists a polar wobble as a resnonsc to each of the components of the tidal potential. 
The magnitude of the polar displacement depends on twe factors: the amplitude of the tidal 
component and its period (in relation to the Chandler period). The maximum periodic contribu- 
tions are: the component number 055.565 with a period of 6798 days (18.613 years) and 50cm 
of pole displacement, this is the component produced by the motion of the lunar ascending node; 
the annual component 056.554 with 37cm of pole displacement and the semiannual 057.555 
with 32 cm. The tidal ' onponents with daily and semidaily periods yield very small pole displace- 
ments. of the order of 0.01 cm. The combined effect of all the periodic components can yield 
as much as 90cm of pole displacement. 

The tidal components with periods longcf than the Chandler period yield larger pole dis- 
placements as the Chandler period is increased, the opposite is true for those components with 
periods shorter than the Chandler period. The changes produced by the ocean tides in the pro- 
ducts of inertia are periodic and regular therefore they can not be the source of excitation of the 
Chandler wobble. 

A possible future line of development is to consider the changes in all the components of 
the inertia tensor (not just 1^2 nonlinear Liouvilie equations by 

means of numerical methods. 
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APPENDIX 1 


INTEGRAL OF THE PRODUCT OF 3 SURFACE HARMONICS 


APPENDIX 1 


INTEGRAL OF THE PRODUCT OF 3 SURFACE HARMONICS 


The surface integral appearing in Equation (3.15) involves the product of three surface spher- 
ical harmonics. In general, 


/ did f Sg'n,' Slm Sgni d\}> = I 

n ft 


I = (_l)M+m (4^) 


Sgm = If(cos0)e>"»'^ 

(S’ + m' )! (L + M)! (« + m)! 


- 11/2 


(C' - m' )! (L - M)! (« - m)! 


/ S' L S \ I a'- L 
\-m' M m/ W 0 


The symbol ( ‘ ■' ) denotes the 3 - j symbol of Wigner, for their evaluation and properties 

\ki kj kj 

see Rotenberg et al., (1959). 
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ANGULAR VARIABLES IN DOODSON’S TIDAL EXPANSION 
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ANGULAR VARIABLES IN DOODSON’S TIDAL EXPANSION 

The tidal expansion is given in terms of six parameters which constitute the frequencies 
appearing in the arguments on the right hand side of Equation (2.3). In terms of their periods 
they are the following: 

V, ; Lunar day, period of 24h 50.47m, 

V 2 : Moon’s mean longitude, period of 27.321582 days, 

V 3 : Sun’s mean longitude, period of 365.242199 days, 

V 4 : Mean longitude of the lunar perigee, period of 8.847 years, 

V 5 : Mean longitude of the lunar ascending node, period of 18.613 years, 

Vg : Mean longitude of perihelion, period of 20,940 years. 

Melchior (1966) gives the following expressions: 

V 2 = 270°.43659 + 481267°.89057T + 0°.00198T2 + 0°.000002j3, 

V 3 = 279°.69668 + 36000°.76892T + 0°.00030T2, 

V 4 = 334° .32956 + 4069°. 03403 T - 0°.01032T2 - 0°.00001 T^, 

V 5 = 259°. 18328 + 1934°. 14201 T + 0°.00208T2 + 0°.000002T3, 

Vg = 281°.22083 + 1°.71902T + 0° .000451^ + 0°.000003T3, 

Vi = t - V 2 + V 3 , 

t: mean solar time, 

T: time expressed in Julian centuries, 

T = 0 at 1899 December 31 , 12h, Om, Os ephemeris time. 
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